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Abstract — Anomaly detection is very useful for hyperspectral 
detection with no a priori information of spectral signature. 
However, most prevailing anomaly detectors are applied directly 
to the all-bands data without considering the characteristics of 
hyperspectral imagery in local spectrum range and some 
computation problem incurred from high dimensional data. In 
this paper, a Dezert-Smarandache Theory (DSmT) based 
anomaly detector is presented to handle these problems. The all- 
bands data are first partitioned into several lower dimensional 
band-subsets. Then the generalized basic belief assignment of 
DSmT reasoning are constructed according to target signal to 
noise ratio (TNR) of different band subsets and probability 
densities of the detection value from different band subsets. 
Experimental results from the real hyperspectral imagery of 
Operative Modular Imaging Spectrometer I (OMIS-I) show that 
our method outperforms the benchmark RX detector (RXD). 
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I. Introduction 

Hyperspectral imagery records the electromagnetic 
spectrum information in high spectral resolution which can be 
utilized to detect target. In hyperspectral anomaly detection, 
we have no a priori information about the spectral signature of 
the target and the background. The most reasonable method is 
to treat pixel whose spectral signature is significantly different 
from that of background as anomaly (target) [1], The most 
commonly used anomaly detector in hyperspectral imagery are 
RX detector (RXD) based on generalized likelihood ratio test 
(GLRT), low probability detection (LPD) and unified target 
detection (UTD) [2-4], etc. All these anomaly detectors are 
based on covariance matrix of all spectral bands (CMASB). 
However, CMASB based anomaly detector has some inherent 
weakness, such as concealing subtle differences between 
materials in some local spectrum ranges, low precision of 
estimating the covariance matrix, and even causing the 
singularity problem in calculation. 

In order to handle these problems, we present a DSmT 
based anomaly detection method in this paper. The whole 
hyperspectral data are first partitioned into several band- 
subsets. This partition reveals the local spectral variation in 
hyperspectral data [5-6]. Also, since the spectral dimension of 
each band-subset is much lower, the singularity problem may 
be avoided when applying CMASB based anomaly detectors 
to each band-subset, and a better estimation of covariance 



matrix can be obtained. After band-subset partition, we apply 
the RXD to each band-subset. Then, the thing remained is to 
determine how to effectively merge the primary detection 
results from individual band-subsets to improve the final 
detection performance. In this paper, we present a novel 
DSmT based anomaly detector to integrate the 
discriminabilities from different band-subsets. 



II. GLRT Based RXD 



Target detection can be stated as a problem of hypothesis 
testing [1][2], The likelihood ratio test is most commonly used 
method for hypothesis testing. However, in some practical 
hyperspectral detection, two likelihood functions of the target 
and the background may depend on some unknown 
parameters. Thus, the generalized likelihood ratio test (GLRT) 
can be utilized as follows 
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where 0 () and 0, are the parameters of the likelihood 



functions of background and target, respectively. x n is the 



observed spectrum of the /7th pixel in hyperspectral imagery, 
and the null hypothesis H 0 and the alternative hypothesis // 
denote target absence and target presence, respectively. Let 
S = [s,,s 2 ,..., Sjy] denote the shape pattern of target and is 



known, t = be the spectrum of target and 

b =[b l ,b 2 ,...,b B ] be the spectrum of background. An 
observed hyperspectral can be expressed as: 



H o :x„ =b + v 
H, :x„ = b + t-s„ +v 
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If b and t follows (V(n 4 ,£ 4 ) and iV(fi,, £,) 
respectively, where covariance matrix of background £ 4 and 
that of target £ ( are identical to £ , and mean vector of 
background \i b is know, we have the hypothesis testing 
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which is multipixel target version of RXD [2][9]. For single- 
pixel target version, (3) is simplified to 
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If E b and (ij are estimated within a neighborhood of 
every test pixel, RXDs shown in (3) and (4) are generalized to 
detectors that are sensitive to local anomaly. 

III. Bands Partition and SNR of Band subsets 

RXD is an anomaly detector directly applied to all-band 
hyperspectral data without considering electromagnetic 
characteristics in local spectrum, sensor noise level and the 
transfer characteristics in imaging. Partitioning the whole 
hyperspectral data into several band-subsets can solve these 
problems, as well as some computational problems. A widely 
used method to partition the complete set of bands into several 
band-subsets is to use correlation coefficients between bands 
[5], because correlation coefficients measure the similarity of 
different bands. For the hyperspectral imagery with B bands, 
the correlation coefficients matrix is defined as 

R oir (i,j) = C(i,j)/JC(i,i)C(j,j) (5) 

for the /th and /th band, where C is the LxL covariance 
matrix with respect to the data matrix Y . Fig. 1 shows the 
correlation coefficients corresponding to all pairs of bands of 
OMIS-1 hyperspectral imagery. The brighter the pixel is, the 
larger the correlation coefficient value is. 

In terms of Fig. 1, the complete set of bands can be 
partitioned into several non-overlapped band-subsets. Since 
correlation coefficients between different bands in same band 
subset are roughly identical, TNRs of bands in each 
band-subset can be viewed as being identical. Assume that 
spectrums of target and background are linearly independent 
and the number of category of these spectrums is less than that 
of spectral band. For covariance matrix of hyperspectral data 
2T in /th band subset, there is[7] 

E.=U,.A,.u;+afl (6) 

where I is the identity matrix, U, is the orthonormal matrix, 
A, is the diagonal and singular matrix and a] is the noise 
level of the ;ith band data. Then we can obtain the TNR 
corresponding to the /th band-subset [7]. 




Figure 1. Correlation coefficient matrix for 80 
bands OMIS-I imagery 



IV. DSMT BASED RX DETECTOR 

The Dezert-Smarandache Theory (DSmT) of plausible and 
paradoxical reasoning proposed in recent years can be 
considered as an extension of the classical Dempster-Shafer 
theory (DST) [8]. Both the probability and the DST are under 
the assumption of exclusivity and exhaustivity of hypotheses. 
Assume there is a discernment frame © = {#,,#,} which is 
composed of two exhaustive elements. Probability deals with 
basic probability assignments (BPA) »/(•) e [0,1] such that 

m(9 l ) + m(9 2 ) = 1 (7) 

The DST, on the other hand, deals with basic belief 
assignments (BBA) »/(•) e [0,1] such that 

m(9 i ) + m(9 2 ) + m(9 l u# 2 ) = l (8) 

The DSmT theory is only under the assumption of 
exhaustivity of hypotheses, and deals with the generalized 
BBA (GBBA) /«(•) e [0,1] such that 

m(9 l ) + m(9 2 ) + m(9 l u 9 2 ) + m(9 l n 9 2 ) = 1 (9) 

One of the distinguished advantages of DSmT is the notion 
of hyper-power set [8]. If let Q = {9 x ,9 2 ,...,9 n } be a finite set 

of n exhaustive elements, hyper-power set D" in the DSmT 
framework is defined as the set of all composite propositions 
built from elements of 0 with union operator u and 
intersection operatorn such that: (1) <p,9 [ ,9 2 ,...,9 n e D e ;(2) 
If A, Be D 0 , then AnBeD 0 and Au Be D e ; (3) No 
other elements belong to D 0 , except those obtained by using 
rule (1) or (2). From 0 , we let the mapping m( ) : D e —> [0,1] 
be as 

m(0) = 0 and E m(A) = 1 (10) 

A<=D e 

The generalized belief and plausibility function are 

Bel(A) = E m(B ) (11) 

B<zA,BeD & 

PI(A)= Z e MB) (12) 

BnA=j>,BeD 0 

Two independent sources of evidences over the same frame 
©with belief functions Bel x (-) and Bel 2 () associated with 
GBBA w,(-) and m 2 () correspond to the conjunctive 
consensus of the sources. It is given by 

VC e D e , m(C) = E mi (A)m 2 (B) (13) 

AnB=C,A,B<=D e 

Each band-subset split from all-bands hyperspectral data 
can be treated as an information source. The data of such an 
information source contains information of local spectral 
signature. Then, the features for target detection can be 
extracted from each source. These features reflect the local 
discriminability of spectral signatures for our detection task. 

Let 9 t and 9 h denote target and background, respectively. 
The value of RXD calculation (before thresholding) of each 
band subset can be viewed as an evidence over the basic frame 
1 9, . 9 b | . Then the corresponding hyper-power set is 
1 9, . 9 b , 9, n 9 b . 9 u 9 h j . Considering that the performance of 
hyperspectral detector is often significantly affected by the 
TNR, we define the GBBAs of target, background, uncertainty 




and paradox corresponding to the zth band subset respectively 
as 



fRXDj 

m i (0 l ) = TNR r j p Rm (x)dx 


(14) 


pRXD: 
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m,{9 pam ) = 0 


(17) 



where 9 uncer is 0, U0 2 , and 9 para is 0, n0 2 , is the 

probability density function (PDF) of RXD value of zth band- 
subset and RXD ( is the RXD value of zth band subset. 
m,(0,) , m,.(0J , m,(9 mcer ) and m,(9 para ) denote supports for 
target, background, uncertainty of target and background, as 
well as paradox of target and background from zth band- 
subset, respectively. Then the support for target from zth and 
y'th band-subset is 



Ri.j = {m i ,m j ) 

= m, (0, )m j (0, ) + m, (0, )m j (9 uncer ) + 1 n ; (0 uncer )m . (0, ) 

+ n h (3 ) m i ( (0 pam )ntj (0, ) + m i (9 uncer )m, ( 0 ;)ora ) 

+m (9 )rn (0 ) 

= w,(0, )m,(0, ) + w,(0, )m . (0„„^) + m i (0 imc Jm j (0 l ) (18) 



where ( , ) denotes the operator corresponding to calculation 
in (13). 

For K ( K >2) band subsets, we utilize the generalized 
fusion rule to construct decision function as follows 

\> A target 

R = 3 , , , ( 19 ) 

[< A background 

to make decision for target detection, where A is the threshold 
for decision. Since {9 t ,9 b ,9 l r\9 b ,9 t U0 fc } is closed under u 



and n operators of {9 r 9 h } , this combination rule of detection 
guarantees that R is a proper generalized belief assignment. 
This combination rale of detection is commutative and 
associative, and can be used for the fusion of detection value 
of K band subset. 



V. EXPERIMENTAL RESULTS 

The data used in this experiment were recorded by OMIS-I 
over Zaoyuan, China in 2001, from the flight height of about 
1,200 meters. The radiance imagery has been calibrated to the 
reflectance imagery. A portion of the image data, i.e., a scene 
of size 100 X 100 pixels is used in our experiments. There are 
several kinds of natural materials and two car target in the 
scene. We name these two targets as target 1 and target 2, 
respectively. 

According to Fig. 1, we partition the all-band data into four 
band subsets: data from the 1st to the 23rd band, data from the 
24th to the 66th band, data from the 67th to the 72ed band, and 
data from the 73rd to the 80th band. These four subband 
datasets are referred as 1st to 4th band-subset respectively. 
Next, we adopt covariance matrix of four band-subsets to 



obtain the TNRs corresponding to these band-subset. These 
TNRs are 155.2974, 158.1240, 223.9299 and 102.2755, 
respectively. 




(c) (d) 

Figure 2. Output of RXD of four band subsets, (a) 1st band subset; (b ) 
2nd band subset; (c) 3rd band subset; (d) 4th band subset. 



The RXD is then applied to each band-subset. Fig. 2(a)- 
2(d) show the RXD results (before thresholding) of the 4 band 
subsets respectively. The sparseness of small targets in the 
whole scene allows the use of the whole hyperspectral data to 
estimate the statistical characteristic the background. 




Figure 3. Detection results of RXD applied to all-band 
hyperspectral data, (a) target 1; (b) target 2. 




Figure 4. Detection results of DSmt based fusion (a) 
target 1; (b) target 2. 











Fig. 3(a)- Fig. 3(b) show detection results of RXD applied 
to full-band data. Fig. 4(a)- Fig. 4(b) show detection results 
utilizing (19). For Fig. 3(a) and Fig. 4(a), the thresholds are set 
by the value that can detect target 1 with the fewest false 
alarms, respectively. For Fig. 3(b) and Fig. 4(b), the thresholds 
are set by the value that can detect target 2 with the fewest 
false alarms, respectively. Compare Fig 3 with Fig. 4, it can be 
clearly seen that our method can detect target 1 and target 2 
with significantly fewer false alarm than the benchmark RXD 
which is applied on all-band data. This demonstrates the 
feasibility and effectiveness of our proposed method. 

VI. Conclusion 

In this paper, we have exploited the DSmT based 
hyperspectral anomaly detection on real OMIS-I data. 
Theoretical analysis and the experimental results have shown 
the effectiveness of our DSmT based detection method and its 
better performance compared with benchmark RXD applied to 
all-band hyperspectral data. 
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